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A theoretical study of the magnetization curves of quasiperiodic magnetic multilayers is presented. 
We consider structures composed by ferromagnetic films (Fe) with interfilm exchange coupling pro- 
vided by intervening nonferromagnetic layers (Cr). The theory is based on a realistic phenomeno- 
logical model, which includes the following contributions to the free magnetic energy: Zeeman, 
cubic anisotropy, bilinear and biquadratic exchange energies. The experimental parameters used 
here are based on experimental data recently reported, which contain sufficiently strong biquadratic 
exchange coupling. 
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I. INTRODUCTION 

The study of the properties of magnetic multilayers 
has been one of the most investigated fields in the last 
decade. The understanding of a number of new and in- 
triguing results became an exciting challenge from both 
theoretical and experimental point of view. In their pi- 
oneer work, Griinberg and colaborators [1] reported evi- 
dences of an antiferromagnetic bilinear exchange coupling 
in Fe/Cr/Fe structures. After that, Baibich et al [2] no- 
ticed a sudden fall in the electrical resistance of Fe/Cr 
magnetic multilayers when an external magnetic field was 
applied. The effect was so striking that was called giant 
magnetoresistance, and recently it has been widely con- 
sidered for applications in information storage technology 
[3]. Through magnetoresistance measurements, Parkin, 
More and Roche [4] observed an oscillatory behavior of 
the exchange coupling in magnetic metallic multilayers 
as a function of the nonmagnetic spacer thickness. This 
work was seminal to a number of experimental studies on 
Fe/Cr/Fe structures with different nonmagnetic spacer 
thickness. Later on, in 1991, Riihrig et al [5] showed 
evidences of a non-colinear alignment (90°) between fer- 
romagnetic layers in Fe/Cr magnetic multilayers, for non- 
magnetic spacer thickness, where the bilinear exchange 
coupling was small. This behavior could not be explained 
considering only the usual bilinear exchange coupling in 
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the free magnetic energy. In fact, the inclusion of a bi- 
quadratic exchange term in the free magnetic energy of 
the system allows the stabilization of non-colinear align- 
ments. Until recently, it was found that the biquadratic 
exchange coupling was too small when compared to the 
bilinear exchange coupling. However, Azevedo et al [6-8] 
presented a number of experimental results in Fe/Cr/Fe 
samples which show the biquadratic exchange coupling 
comparable to the bilinear exchange coupling. There- 
fore, the biquadratic coupling can play a remarkable role 
in the properties of magnetic multilayers. 

On the other hand, from an experimental point of view, 
due to the rapid development of the crystal growth tech- 
niques, is now possible to tailor a wide class of magnetic 
multilayers, whose film thickness is extremely well con- 
trolled. As a consequence, there are new magnetic phases 
and properties which are not shared by the constituent 
materials. 

It is known that the magnetic properties can depend 
strongly on the stacking pattern of the layers. Under this 
aspect, the physical properties of a new class of artificial 
material, the so-called quasiperiodic structures, became 
recently an attractive field of research. Quasiperiodic 
structures, which can be idealized as the experimental re- 
alization of a one-dimensional quasicrystal, are composed 
by the superposition of two (or more) building blocks that 
are arranged in a desired manner. They can be defined as 
an intermediate state between an ordered system (a peri- 
odic crystal) and a disordered one (an amorphous solid) 
[9,10]. One of the most interesting features of these sys- 
tems is that the long range correlations, induced by the 
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construction of the systems, are reflected in their various 
spectra. In fact, many physical properties of quasiperi- 
odic systems have been studied such as light propagation 
[11], phonons [12], electronic transmission [13], polari- 
tons [14], and magnons [15] . In all of these situations, 
despite the diversity of the systems, a common feature 
is present, namely, a fractal spectra of energy, which can 
be considered as their basic signature [16,17]. However 
only very recently some efforts were done towards the un- 
derstanding of the properties of quasiperiodic magnetic 
multilayers [15,18]. 

The main aim of this paper is a contribution to the 
understanding of the effects of the quasiperiodic arrange- 
ment on the magnetization curves in magnetic multilay- 
ers. We are interested in new magnetic phases and align- 
ments that are only due to the quasiperiodicity of the 
system. We have studied Fe/Cr(100) structures which 
follow a Fibonacci and a double period or generalized Fi- 
bonacci quasiperiodic sequences. 

The layout of the paper is as follow: In section II we 
discuss the physical model used here, with emphasis in 
the description of the quasiperiodic sequences. In section 
III we define the contributions to the magnetic energy. 
The numerical methods, used to obtain the equilibrium 
configuration, are described in section IV. In section V, 
the results are presented and discussed. Finally, we draw 
the conclusions in section VI. 



II. PHYSICAL MODEL 

A quasiperiodic structure can be experimentally con- 
structed juxtaposing two building blocks (or, as consid- 
ered here, building layers) following a given quasiperiodic 
sequence. We choose Fe as the building layer associated 
with the letter A, and Cr as the building layer associated 
with the letter B (see Fig. 1). Therefore, we only take 
into account generation of sequences that start and finish 
with an Fe building layer, which means an even number 
of Fe layers, to guarantee a real magnetic counterpart. 
In this way we avoid also the intriguing behavior found 
when even and odd numbers of Fe layers are considered, 
as discussed in [19]. In this paper we have considered two 
quasiperiodic sequences, namely, the Fibonacci and the 
double period sequences. 



A. The Fibonacci magnetic multilayers 

The Vth generation of the Fibonacci sequence can be 
determined appending the N — 2 generation to the N — 1 
one, i.e., Sn = SW-1SW-2 (V > 2). This algorithm 
construction requires initial conditions which are chosen 
to be So = B and Si = A. The Fibonacci generations 
can also be alternatively obtained by an iterative process 
from the substitution rules (or inflation rules), A — » AB, 
B — > A. The Fibonacci generations are: 



So = [B] , Si = [A] , S 2 = [AB] , S 3 - [ABA] etc. 

In a given generation Sn, the total number of letters is 
given by the Fibonacci number Fn, which is obtained by 
the relation F N — F N -i + F N -2, with F a — Fi = 1. 
Also, Fn-i and Fn-2 are the number of letters A and 
B, respectively. As the generation order increases (N » 
1), the ratio Fn/Fn-i approach to r = (1 + \/5)/2, an 
irrational number which is known as the golden mean. It 
is also possible to obtain the number of letters A and B 
for a given generation by the substitution matrix of the 
Fibonacci sequence Mp from [20], 
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Here (n A + ,n B + ) are the number of letters A and B in 
the (N+l)th generation, and (n A , rig) are the number of 
letters A and B in the Nth generation. The explicit form 
of the substitution matrix for the Fibonacci sequence is, 



M F = 



1 1 
1 



(2) 



whose first eigenvalue, A, is the golden mean r. 

In Fig. 1, we show the third and fifth Fibonacci gen- 
erations and their magnetic counterparts. Note that 
the third Fibonacci generation corresponds to a trilayer 
Fe/Cr/Fe, and in the fifth Fibonacci generation there is 
a double Fe layer. It is easy to show that the Fibonacci 
magnetic multilayers, for any generation, are composed 
by single Cr layers, single Fe layers and double Fe layers. 
The number of Fe single layers is 1 + Fn-2, the number 
of Fe double layers is — 1 + Fn- 1 — -Fjv-2 and the number 
of Cr layers is F N -2- It should be observed that only 
odd Fibonacci generations have a magnetic counterpart 
(they start and finish with an Fe building layer). 



B. The double period magnetic multilayers 

The Nth generation of the double period sequence can 
be obtained from the relations, 



with 



Sat — Sjv-iS^-l, 



S N = S N . 1 S N -i (V>2), 



(3) 



(4) 



The initial conditions are So = A e Si = AB. We can 
alternatively use the substitution rules A — > AB, B — > 
AA. The double period generations are: 

S = [A] , Si = [AB] , S 2 = [ABAA] , etc. 
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In a given generation Sn, the total number of letters is 
2^, and the number of letters A and B for consecutive 
generations can be related by the substitution matrix of 
the double period sequence Md p , i.e.: [20] 



E T = E z + E ca + Eu + E bq 



(7) 
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are the number of letters A and B in 



the (N + l)th generation, and (n 



N N\ 
A > n B , 



are the number 



of letters A and B in the Nth generation. As the gener- 
ation number increases (N » 1), the ratio between the 
number of letters A and B tends to 2. The explicit form 
of the substitution matrix for the double period sequence 
is, 



M F = 



1 1 

2 1 



(6) 



In Fig. 2 we show the second and fourth double period 
generations and their magnetic counterparts. The dou- 
ble period magnetic multilayers are composed by single 
Fe layers, double Fe layers, triple Fe layers and single Cr 
layers. It should be observed that, contrary to the Fi- 
bonacci case, only even double period generations have a 
magnetic counterpart. 



III. MAGNETIC ENERGY 

We consider magnetic multilayers whose constituents 
are Fe ferromagnetic films, separated by Cr non-magnetic 
films. We take the xj/-plane as the film plane and the z- 
axis as the growth direction. We consider that the mag- 
netic films are uniformly magnetized and that they be- 
have as monodomains. We also consider that they do not 
present dynamical excitations and that the very strong 
demagnetization field, generates by tipping the magneti- 
zation out of the plane, will suppress any tendency for 
the magnetization to tilt out of plane. Therefore, the de- 
grees of freedom of the magnetizations are restricted to 
the aiy-plane. The interfilm exchange couplings between 
the ferromagnetic films are weak when compared to the 
strong exchange couplings between spins within a given 
ferromagnetic film. Therefore, we can represent the ferro- 
magnetic films as classical magnetizations M, composed 
by the real spins within the films, which are strongly cou- 
pled by the intrafilm exchange coupling. These classical 
magnetizations interact through the interfilm exchange 
coupling and they can present some anisotropy depend- 
ing on the structure studied. It should be noted that 
this system is isomorphous to a one-dimensional chain of 
classical spins. 

The global behavior of this system is well described by 
a realistic phenomenological theory in terms of the free 
magnetic energy [7], i.e., 



Here E z is the Zeeman energy (between the ferromag- 
netic films and the external applied magnetic field), E ca 
is the cubic crystalline anisotropy energy (which we con- 
sider present in the ferromagnetic films) and Eu and Ey, q 
are the bilinear and the biquadratic exchange coupling 
energies (between the ferromagnetic films), respectively. 

The explicit form of the free magnetic energy can be 
written as, 

n 

E T = -Y^UMi-E 



f K 

i—i I l \ 



^ Mi- M l+1 \X (Mi ■ M t+1 f 



\Mi\\M, 



t+ll 



\Mi\>\M i+ tf 



(8) 



Here, H is the external magnetic field which is applied 
in the film plane, U is the thickness of the ith Fe layer, 
Mj is the classical magnetization of the ith Fe layer, and 
K ca is the cubic anisotropy constant. Also, Ju and Jb q 
are the bilinear and the biquadratic exchange couplings, 
respectively. This expression, after a tedious but straight 
calculation, takes the form, 

n 

e t = Y, {-UMiH oaa(6i - H ) + -t l K ca sm 2 (20 i )} 

i=l 

n-l 

+ 2 {~ Jbl C ° S ( 6i ~ dl+1 "> + J bl COs2 (^ - fc+l)}- ( 9 ) 

i=l 

Here 9i is the angular orientation of the magnetization 
of the ith Fe layer and Oh is the angular orientation of 
the magnetic field. From this point we consider 9h = 0, 
which means that the magnetic field is applied along the 
easy axis. It is usual to write the total free magnetic 
energy in terms of experimental parameters, like 



Hca — 

H bl = 
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(11) 

(12) 



In this way we obtain a final expression for the free mag- 
netic energy per unit area, 



Ej> 
tM~ s 



n 1 

J2(U/t){-H cos(0i) + -H ca sm 2 (20A} 
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+ J2 i-Hbi cos(0 4 - + H^cos 2 ^ - l+1 )}. (13) 
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Here t is the thickness of a single Fe layer which is con- 
sidered to be the basic tile, Mi is assumed to be equal 
to Ms (the saturation magnetization), and H ca is the 
cubic anisotropy field which turns the (100) direction an 
easy direction. Hu is the bilinear exchange coupling field 
which favors antifcrromagnetic alignment when negative, 
and ferromagnetic alignment when positive. Hb q is the 
biquadratic exchange coupling field which is experimen- 
tally found to be positive and favors a non-colinear align- 
ment (90°) between two adjacent magnetizations. 

Once the free magnetic energy is determined, we can 
calculate the equilibrium configuration for specific values 
of the experimental parameters as a function of the ex- 
ternal applied field. In simple situations, the equilibrium 
configuration can be analytically obtained by equating 
to zero the derivatives of the magnetic energy with re- 
spect to the angle 9. However, in most cases this leads to 
transcendental equations which can not be analytically 
solved. From a numerical point of view, many methods 
have been proposed to calculate the equilibrium positions 
of the magnetizations. In the next section we describe the 
methods used in this paper. 



IV. NUMERICAL METHODS 

In this section, we want to find the global minimum of 
the cost function, 



Et — Et (0i , 62 



(14) 



where 9 n can assume values in the range [0, 27r] and it 
defines a n-dimensional space. When the dimension of 
this space is high, the cost function has a rough surface, 
i.e., there are many local minima which make difficult 
to find the global minimum. There are many numerical 
methods to solve this problem [21]. In our specific case 
two methods were successfully used, namely, simulated 
annealing and the so-called gradient method. 



A. Simulated annealing method 

Introduced by S. Kirkpatrick [22], simulated annealing 
(SA) comes from the fact that the heating (annealing) 
and slowly cooling a metal, brings it into a more uni- 
forme crystalline state, which is believed to be the state 
where the free energy of bulk matter takes its global min- 
imum. The role played by the temperature is to allow the 
configurations to reach higher energy states with proba- 
bility given by Boltzmann's exponential law. Then they 
can overcome energy barriers that would otherwise force 
them into local minima. In general, a simulated anneal- 
ing technique can be written as follows: 

(i) Choose an initial point in parameter space, corre- 
sponding to an initial configuration {9}j, and cal- 
culate the associated energy Ej. 



(ii) Choose a second point in parameter space, corre- 
sponding to a second configuration {9}j + i, and cal- 
culate the associated energy Ej + \. 

(iii) If AE = Ej+i - Ej < 0, {9}j+i is the new config- 
uration of the system. 

(iv) If AE > 0, we define the probability p = 
exp(— AE/ksT) and choose a random number < 
x < 1. If x > p, {9}j+i is the new configuration of 
the system. Otherwise, {9}j is maintained as the 
configuration of the system. 

(v) This procedure is executed again and again until 
the equilibrium is reached. 



B. The gradient method 

The second method that we have used was the so-called 
gradient method [23]. This method is based on the di- 
rectional derivative of the cost function (the magnetic 
energy) in the search of its global minimum. In this way, 
we need to calculate the gradient of Et with relation to 
the set {9}, 



VE T = 



dE 7 
~d9~ 
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From this relation we execute the following algorithm to 
find the equilibrium configuration, 

(i) We generate a configuration in the parameter space 
{9}j from which we calculate the associated energy 
Ej and the gradient of the cost function. 

(ii) A second point in the parameter space is generated 
by {9}j+i = {6}j — aVEr- Here a controls the size 
of the displacement in the direction — VEt- 

(iii) The energy of the second point is calculated and if 
Ej + i > Ej, the parameter a (the size of the dis- 
placement) is divided by two and we go back to (ii). 
Otherwise, we instead generate a new configuration 
from {9}j + i. 

In the last step the reduction of a is limited by the 
precision value required for e. This limit is reached when 

aVE T < e. 

We have used the two methods discussed above to 
obtain the equilibrium positions of the magnetizations. 
Each method was applied for each value of the applied 
magnetic field and for each set of experimental parame- 
ters. We choose the configuration with the lowest energy 
furnished by both methods as the equilibrium configura- 
tion. 
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V. NUMERICAL RESULTS 

In this section we present the numerical results ob- 
tained for the magnetization curves of quasipcriodic mag- 
netic multilayers. In all situations we have considered the 
cubic anisotropy effective field H ca = 0.5 kOe which cor- 
responds to Fe(100) with t > 30A. In our calculations 
we have used two sets of experimental values for the bi- 
linear and biquadratic exchange coupling: (i) the first 
one with Hu = — 1.0 kOe and H^q = 0.1 kOe. It lies 
in the region of the first antiferromagnetic peak of the 
bilinear exchange coupling, corresponding to a realistic 
sample whose Cr thickness is about 10A; (ii) the second 
set with H b i = -0.035 kOe and H bq = 0.035 kOe. It is 
in the region of the second antiferromagnetic peak of the 
bilinear exchange coupling, corresponding to a realistic 
sample whose Cr thickness is about 25A. 



A. Fibonacci magnetic multilayers 

The magnetization curves for the first set of parameters 
of the Fibonacci magnetic multilayers are shown in Fig. 3. 
For the third generation (which corresponds to the well 
known Fe/Cr/Fe trilayer), in the low field region, the 
magnetizations are antiparallel. As the field increases, 
they continuously rotate toward the field direction (sec- 
ond order phase transition) and the saturation is reached 
when the external magnetic field H ~ 1.91 kOe. For the 
fifth generation there are two first order phase transitions 
at H ~ 0.71 kOe and H ~ 0.87 kOe, respectively. The 
saturation is reached at H <~ 2.93 kOe. For the seventh 
generation, there are three first order phase transitions at 
H ~ 0.28 kOe, H ~ 0.96 kOe and H ~ 1.06 kOe, respec- 
tively. The saturation is reached at H ~ 3.03 kOe. For 
this set of parameters the majority of the transitions are 
of second order. Note that due to the different Fe layer 
thickness, for the fifth and seventh generations, the mag- 
netization is not zero even for zero magnetic field. In Fig. 
4 we show the results for the second set of parameters. 
For the third generation, due to the strong biquadratic 
field, there is no antiparallel phase in the low field region. 
Two magnetic phases are present: 90° (H < 72 Oe) and 
saturated (H > 72 Oc). The fifth generation presents 
three magnetic phases: (i) 90° (H < 72 Oe); (ii) almost 
saturated (72 Oe < H < 0.14 kOe); and (iii) saturated 
(H > 0.14 kOe). The seventh generation presents four 
magnetic phases from 90° (H < 36 Oe) to the saturated 
regime (H > 0.14 kOe). All transitions are of first or- 
der. Note the striking self-similar pattern shown by the 
magnetization profile in this figure (see the windows). 

B. Double period magnetic multilayers 

Figs. 5a and 5b show our results for the double period 
magnetic multilayers, considering the first set of param- 



eters. For the second generation, due to the double Fc 
layer, the magnetization has about 1/3 of its saturation 
value for zero magnetic field. There is a first order phase 
transition from antiparallel to an asymmetric phase at 
H <~ 0.69 kOe. In this phase, the magnetizations arc 
asymmetrically oriented along the magnetic field. When 
H ~ 1.34 kOe the saturated phase emerges. For the 
fourth generation, for zero magnetic field, the magneti- 
zation has about 10% of its saturation value due to the 
different thickness of the Fe layers. There is a first or- 
der phase transition at H <~ 0.29 kOe and the saturation 
is reached at H ~ 3.27 kOe. All other phase transitions 
are of second order. For the second set of parameters (see 
Fig. 6), on the contrary, all transitions arc of first order. 
For this set of parameters, there is no antiparalell phase 
in the low field region, due to the strong biquadratic 
field. For the second generation, the magnetization is 
about 2/3 of its saturation value when H = 0. There 
are two magnetic phases: (i) 90° (0 < H < 72 Oc) and 
(ii) saturated (H > 72 Oc). For the fourth generation, 
the magnetization for zero magnetic field is about 1/2 of 
its saturation value. Four magnetic phases are present, 
from the 90° (0 < H < 38 Oe) to the saturated phase 
(H > 0.14 kOe). As in the Fibonacci self-similar 
pattern is also present in the magnetization curves (see 
the window). 

VI. CONCLUSIONS 

We have studied quasiperiodic magnetic multilayers, 
composed by ferromagnetic Fe layers separated by non- 
magnetic Cr layers, arranged according to the Fibonacci 
and double period quasipcriodic sequences. We consider 
that the Fe layers are linked by bilinear and biquadratic 
exchange couplings through Cr layers and present cubic 
anisotropy. The external magnetic field is applied in the 
plane of the layers and along an easy axis. We have used 
two numerical methods to determine the equilibrium con- 
figurations of the layers's magnetizations. The magneti- 
zation curves of these artificial structures were calculated 
considering two sets of experimental parameters recently 
reported [6-8] . Our results show that quasiperiodic mag- 
netic multilayers exhibit a rich variety of configurations 
induced by the external magnetic field. In particular 
two points may be emphasized: (i) the effect of different 
thickness of Fe layers and (ii) the effect of the biquadratic 
exchange coupling. 

The effect of different thickness of Fe layers is evident 
in the low field region. In that region, due to these differ- 
ences, there is a net magnetization even if the alignment 
is antiparallel and the external magnetic field is zero. Be- 
sides, the nature of the phase transitions are changed by 
the different thickness (Fig. 3a shows only second order 
phase transitions, while Fig. 5a presents an additional 
first order phase transition). These results suggest that, 
varying the thickness of Fe layers, it is possible to tai- 
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lor magnetic multilayers to present desired specific phase 
transitions and critical fields. However, as the thickness 
of Fe layers increases, the crystalline anisotropy of Fc 
(100) films on Cr (100) also increases. Fortunately, as 
a characteristic of the quasiperiodic multilayers arrange- 
ments considered here, the maximum number of joint Fe 
layers is two (for the Fibonacci case) and three (for the 
double period case), no matter is the value of their gen- 
eration numbers. Besides, from a thickness greater than 
40A, the crystalline anisotropy reaches saturation [7]. 

On the other hand, the biquadratic exchange cou- 
pling plays a remarkable role in the features of the mag- 
netization curves. For example, when the bilinear ex- 
change coupling prevails, the majority of the transitions 
are of second order character (see Figs. 3 and 5). How- 
ever, when the biquadratic exchange is compared to the 
bilinear one, in the presence of a stronger crystalline 
anisotropy [8], the transitions are characterized by dis- 
continuous jumps in the magnetization that indicate first 
order phase transitions. This can be considered as the 
basic signature of the biquadratic exchange coupling (see 
Figs. 4 and 6), although for the case where there is no 
biquadratic term, a first order phase transition appears 
due to the anisotropy [24]. Furthermore, as shown by 
the windows in these figures, the magnetization curves of 
higher generations reproduce the magnetization curves 
of lower generations. This self-similar behavior is a gen- 
eral characteristic of quasiperiodic systems, although it is 
not present when the bilinear exchange prevails (Figs. 3 
and 5). A possible explanation for these different behav- 
iors is because the biquadratic exchange coupling induces 
long range correlations that emphasize the quasiperiod- 
icity of the system. These long range correlations make 
the whole structure seeing its quasiperiodicity, which is 
reflected in the magnetization curves. This argument is 
reinforced by previous works on the correlation lengths of 
magnetic systems presenting biquadratic exchange cou- 
pling (see, for example, S0rensen and Young [24]). 

The most appropriate experimental technique for 
studying the magnetization curves of magnetic films is 
the magneto-optical Kerr effect (MOKE) [8]. However, 
because the MOKE measurements provide surface sensi- 
tivity on the scale of the optical penetration depth (<~ 10 
A), it is necessary to use also a superconductor quan- 
tum interface device (SQUID) magnetometry [19]. The 
two techniques prove complementary in understanding 
the switching behavior of the multilayer films, as far as 
the magnetization curves are concern. We hope that 
the present results can stimulate experimental studies of 
these structures. 
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Figure Captions 



1 . The third and fifth Fibonacci generations and their 
magnetic counterpart. 

2. Same as Fig. 1 for the second and fourth double 
period generations. 

3. Magnetization versus applied field for the third 
(a), fifth (b) and seventh (c) Fibonacci generations 
with \Hb q \/\Hu\ = 0.10, corresponding to a real- 
istic sample whose Cr thickness is about 10A. We 
have considered the cubic anisotropy effective field 
H ca = 0.5 kOe, which corresponds to Fe(100) with 
t > 30A. 

4. Same as Fig. 3, but for \Hb q \/\Hbi\ = 1.0, corre- 
sponding to a realistic sample whose Cr thickness 
is about 25A. 

5. Magnetization versus applied field for the second 
(a) and fourth (b) double period generations with 
\Hb q \/\Hbi\ = 0.10. The cubic anisotropy effective 
field is again H ca = 0.5 kOe. The Cr thickness is 

10A. 

6. Same as Fig. 5, but for \Hb q \/\Hbi\ — 1.0, corre- 
sponding to a realistic sample whose Cr thickness 
is about 25A. 
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